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We study numerically the SU(2) Landau gauge transverse and longitudinal gluon propagators at 
non-zero temperatures T both in confinement and deconfinement phases. The special attention is 
paid to the Gribov copy effects in the IR-region. Applying powerful gauge fixing algorithm we find 
that the Gribov copy effects for the transverse propagator Dt(p) are very strong in the infrared, 
while the longitudinal propagator Dl(p) shows very weak (if any) Gribov copy dependence. The 
value Dt{0) tends to decrease with growing lattice size; however, Dt{0) is non-zero in the infinite 
volume limit, in disagreement with the suggestion made in We show that in the infrared region 
Dt{p) is not consistent with the pole- type formula not only in the deconfinement phase but also 
for T < T c . We introduce new definition of the magnetic infrared mass scale ('magnetic screening 
mass') hum- The electric mass vtie has been determined from the momentum space longitudinal 
gluon propagator. We study also the (finite) volume and temperature dependence of the propagators 
as well as discretization errors. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38. Aw 
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I. INTRODUCTION 



One of the most interesting features of the quantum 
chromodynamics at finite temperature is the transi- 
tion from the confinement to the deconfinement phase. 
This transition separates a low-temperature phase, 
which is expected to be highly non-perturbative and 
characterized by quark and gluon confinement, from 
a high-temperature - quark-gluon plasma (QGP) - 
phase, where color charges should be deconfined. The 
conjecture of the existence of the QGP has been sup- 
ported by recent observations of the collective effects 
in ultrarelativistic heavy-ion collisions at SPS and 
RHIC (see, e.g., Q and references therein). 

The non-perturbative - first principle - calculation 
of gauge variant gluon (as well as ghost) propagators 
is of interest for various reasons. These propagators 
are expected to show different behavior in each phase 
and, therefore, to serve as a useful 'order parameters', 
detecting the phase transition point T c . One expects 
that their study can shed the light on the mechanism 
of the confinement-deconfinement transition. Another 
reason is that for the reliable phenomenological anal- 
ysis of high-energy heavy-ion collision data, it is im- 
portant to obtain information on the momentum de- 
pendence of the longitudinal (electric) gluon propa- 
gator Dl^p) and transverse (magnetic) gluon propa- 
gator Dt(j>), especially in the (deep) infrared region. 
One example is the study of the radiative energy loss 
in dense nuclear matter (jet quenching) which results 
from the energy loss of high energy partons moving 



through the plasma (see, e.g., [MZj). Also, the non- 
perturbatively calculated lattice propagators are to 
be used to check the correctness of various analyti- 
cal methods in QCD, e.g., Dyson-Schwinger equations 
(DSE) method. For study of the gluon propagator us- 
^ing DSE approach at finite temperature see e.g. 0- 

The lattice study of the finite temperature SU(2) 
gluon propagators in Landau gauge has been 
formed in a number of papers (see, e.g., papers 
[lH). In Ref. [H| the electric and magnetic propaga- 
tors were studied in both coordinate and momentum 
spaces and in 4 and 3 dimensions. The conclusion 
was made that the magnetic propagator had a com- 
plicated infrared behavior which was not compatible 
with simple pole mass behavior. It was also found 
that this propagator had strong volume and gauge 
dependence. In Refs. [13, EH results for both gluon 
and ghost propagators in momentum space were pre- 
sented, but some important questions, e.g. Gribov 
copies effects, infrared behavior and scaling were not 
addressed. 

In paper [l| it has been suggested that the proxim- 
ity of the Gribov horizon at finite temperature forces 
the transverse gluon propagator Dt(p,Pa = 0) to 
vanish at zero three-momentum. If this is the case, 
then the finite-temperature analog of Gribov formula 
\p\ 2 /(\p\ 4 + Mj,) = l/{\p\ 2 +ml ff (p)) suggests that 
the effective magnetic screening mass m e //(p*) be- 
comes infinite in the infrared (interpreted as a mag- 
netic gluons 'confinement'). 

The Gribov copy effects still remain one of the most 
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serious problem in the lattice calculations, at least, in 
the deep IR-region. In our study we employ the gauge 
condition which requires the Landau gauge fixing 
functional F (sec the definition in Section [TTJ) to take 
extrema as close as possible to the global extremum. 
This choice for the gauge condition is supported by 
the following facts: a) a consistent non-perturbative 
gauge fixing procedure proposed by Parrinello-Jona- 
Lasinio and Zwanziger (P JLZ-approach) [I?} pre- 
sumes that the choice of a unique representative of the 
gauge orbit should be through the global extremum 
of the chosen gauge fixing functional; b) in the 
case of pure gauge U(l) theory in the weak coupling 
(Coulomb) phase some of the gauge copies produce 
a photon propagator with a decay behavior inconsis- 
tent with the expected zero mass behavior (l8l - [2(l |. 
The choice of the global extremum permits to obtain 
the physical - massless - photon propagator. 



In a series of papers [21H25I ] we investigated the Gri- 
bov copy effects in the Landau gauge gluon (and/or 
ghost) propagators in zero temperature SU(2) gluo- 
dynamics. It has been demonstrated unambiguously 
that these effects are strong in the infrared 1 . Thus 
the Gribov copy effects reduction is very important for 
the infrared behavior studies. Recently it has been 
pointed out in Ref. (2(| that lattice results for the 
infrared gluon and ghost propagators free of Gribov 
copies effects would help to discriminate between scal- 
ing and decoupling solutions of the Dyson- Schwinger 
equations. In this paper we undertake a careful study 
of the Gribov copy effects in the Landau gauge gluon 
propagator at finite temperature. We employ the 
gauge fixing procedure which we used recently in our 
study at zero temperature [25[ with changes dictated 
by nonzero temperature (see Section [TTT]) . 



II. GLUON PROPAGATORS: THE 
DEFINITIONS 

We study the SU(2) lattice gauge theory with the 
standard Wilson action 



where f3 = 4/<?q and go is a bare coupling constant. 
The link variables U xfl € SU(2) transform under 
gauge transformations g x as follows: 



U x » A E7*. - gtU x ,g x+f , ; g x e SU(2) . (1) 

Our calculations were performed on the asymmetric 
lattices with lattice volume V — L4 ■ L i sl where L4 is 
the number of sites in the Ath direction. The temper- 
ature T is given by 



T = 



(2) 



where a is the lattice spacing. We employ the stan- 
dard definition of the lattice gauge vector potential 
Ac+H/2,n [13]: 

The Landau gauge fixing condition is 



Ac+/t/2,/i n ■ (3) 



{dA) x = J2 (A+/i/2; j u - = , (4) 



Also we attempt to make a careful analysis of (fi- 
nite) volume and temperature dependence of Dt and 
Dl as well as of scaling violations. 



which is equivalent to finding an extremum of the 
gauge functional 



W ^ 2 

Xfl 



(5) 



Section |TT] contains main definitions as well as some 
details of simulations and gauge fixing procedure wc 
use. Section HIIl is dedicated to the study of the Gribov 
copy effects. Volume and temperature dependence of 
the propagators as well as discretization errors are dis- 
cussed in Section [TV] Section [V] is dedicated to the 
discussion of the screening masses and Section IVII is 
reserved for conclusions and discussion. 



with respect to gauge transformations g x . After re- 
placing U => U 9 at the extremum the gauge condition 
((4]) is satisfied. 

The (unrenormalized) gluon propagator D"„(p) is 
defined as follows 



D a Jl(p) = -(A;(k)Al(-k)) 
9o 



1 Unfortunately, authors of [15|| cite our paper 22] in a com- 
pletely misleading context. 



where A(k) represents the Fourier transform of 
the gauge potentials defined in Eq.([3|) after hav- 
ing fixed the gauge, e (—L s /2, L s /2] and k 4 6 
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(—£4/2,-^4/2]. The physical momenta p^ are given 
by api — 2 sin (nki/L s ), ap 4 = 2 sin (irk^/L^). 

In what follows we consider only soft modes p 4 = 0. 
The hard modes (j>4 7^ 0) have an effective thermal 
mass 2irT\k4\ and behave like massive particles 2 . 

As is well known, on the asymmetric lattice there 
are two tensor structures for the gluon propagator 

M ■■ 

DfM = S ab (pT v {p)D T {p) + P^(p)D L (p)) , (6) 

where (symmetric) orthogonal projectors P^ L (p) are 
defined at p = (p ^ 0; P4 = 0) as follows 

p iM = (V^|r)> P »= ! (7) 

P 4 L 4 (p) = 1 ; i$(p)=0. (8) 

Therefore, two scalar propagators - longitudinal 
D l (p) and transverse Dt {p) - are given by 



d t (p) = gEE^W: 

a— 1 2—1 

1 3 

6 a=l 

(9) 

For ;? = propagators -Dt(O) and -Dl(O) are defined 
as follows 

3 3 

a—1 i—l 

D l(0) = ^E^o(O). (10) 

6 a=l 

The transverse propagator Dt(p) is associated to 
magnetic sector, and the longitudinal one Dl(jp) - to 
electric sector. 

We generated ensembles of up to two thousand inde- 
pendent Monte Carlo lattice field configurations. Con- 
secutive configurations (considered as independent) 
were separated by 100 (for L s < 32) or 200 (for 
L s > 32) sweeps, each sweep being of one local heat- 
bath update followed by L s /2 microcanonical updates. 
In Table U we provide the full information about the 
field ensembles used throughout this paper. 



2 Let us note that the 4th euclidian component p4 7^ has no 
physical meaning. 



In order to keep finite- volume effects under control, 
we considered a few different lattice volumes for each 
temperature. The choice of the 6 x 48 3 lattice at = 
2.635 is important for the check of the scaling behavior 
(see Section HV|) . 
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800 
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3.596 


0.055 


6 


40 


2.0 


800 


40 


2.635 


3.596 


0.055 


6 


48 


2.0 


800 


40 



TABLE I: Values of /3, lattice sizes, temperatures, num- 
ber of measurements and number of gauge copies used 
throughout this paper. To fix the scale we take ^/a = 440 
MeV. 

For gauge fixing we em ploy the Z(2) flip operation 
as has been proposed in [231 ]. It consists in flipping 
all link variables U xfl attached and orthogonal to a 3d 
plane by multiplying them with — 1. 

Such global flips are equivalent to non-periodic 
gauge transformations and represent an exact symme- 
try of the pure gauge action. The Polyakov loops in 
the direction of the chosen links and averaged over the 
3d plane obviously change their sign. At finite temper- 
ature we apply flips only to directions fi — 1,2, 3. In 
the deconfincmcnt phase, where the Z{2) symmetry 
is restored, the Z(2) sector of the Polyakov loop in 
the [i = 4 direction has to be chosen since on large 
enough volumes all lattice configurations belong to 
the same sector, i.e. there are no flips between sec- 
tors. We choose sector with positive Polyakov loop. 
In the confinement phase one may use a flip in the 
fj, = 4 direction. However, in a test run we have found 
that at /3 = 2.26 studied in this paper the maximal 
gauge fixing functional ([S]) has been found in the pos- 
itive Polyakov loop sector in more than 90 % cases. 
To save computer time we stick to this sector for all 
configurations at this f3. Therefore, in our study the 
flip operations combine for each lattice field configu- 
ration the 2 3 distinct gauge orbits (or Polyakov loop 
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sectors) of strictly periodic gauge transformations into 
one larger gauge orbit. 

Following Ref. 25] in what follows we call the com- 
bined algorithm employing simulated annealing (SA) 
algorithm (with finalizing overrelaxation) and Z(2) 
flips the 'FSA' algorithm. For every configuration 
the Landau gauge was fixed N copy = 24(40) times 
(3(5) gauge copies for every flip-sector) on lattices 
with L s < 26(L S > 32), each time starting from a 
random gauge transformation of the mother config- 
uration, obtaining in this way N copy Landau-gauge 
fixed copies. We take the copy with maximal value of 
the functional (O as our best estimator of the global 
maximum and denote it as best ("&c ") copy. In or- 
der to demonstrate the Gribov copy effect we compare 
with the results obtained from the randomly chosen 
first ("/c ") copy. 

We present results for the unrenormalized propa- 
gators and make comments on renormalization when 
propagators computed at different values of /3 are com- 
pared or comparison with the renormalized results of 
other groups is necessary. Other details of our gauge 
fixing procedure are described in our recent papers 

To suppress 'geometrical' lattice artifacts, we have 
applied the "a-cut" [29|, i.e. Trki/L s < a , for every 
component, in order to keep close to a linear behav- 
ior of the lattice momenta pi ~ (2nki)/(aL s ), ki G 
(— L s /2, L s /2\. We have chosen a = 0.5. Obviously, 
this cut influences large momenta only. We did not 
employ the cylinder cut in this work. 



III. GRIBOV COPY EFFECTS AND LARGE 
L s BEHAVIOR OF D T (0) 

As has been already pointed above, the Gribov copy 
problem still remains acute, at least, in the deep in- 
frared region and the choice of the efficient gauge fix- 
ing method is very important. The importance of this 
choice is demonstrated in Fig.[T]taken from our recent 
paper [25[ . In this Figure we compare our be FSA re- 
sults for the bare gluon propagator D(p) calculated on 
a 44 4 lattice with those of the standard fc OR method 
obtained for an 80 4 lattice and also with the fc SA re- 
sults. In particular, we observe that the OR method 
with one gauge copy produces completely unreliable 
results for the range of momenta |p|^$Q.7 GeV. (The 
detailed discussion can be found in [25j]) 

Let us define the normalized difference of the fc and 
be transverse propagators At(p) : 



A T (p) 



D f T c (p) ~ D b T °(p) 



D b T c {p) 



(11) 
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"3 6 
Q 
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O 
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• 


44 4 ; 
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where the numerator has been obtained by averaging 
over all configurations of the difference between fc and 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Ipl; Gev 

FIG. 1: Comparison of data obtained for be FSA gauge fix- 
ing with those obtained with the standard fc OR - method 
and the fc SA algorithm (this Figure from paper [25[) 



be transverse propagators calculated for every config- 
uration, this average being normalized to the be (av- 
eraged) transverse propagator. In a similar way one 
can define also the normalized difference of the fc and 
be longitudinal propagators A^(p). 

The longitudinal propagators Dl(p) demonstrate 
very weak dependence on the choice of Gribov copy. 
In Figure [5] we show the momentum dependence of 
A L (p) on the 4 x 48 3 lattice at (3 = 2.35, 2.512 and 
f3 = 2.26, i.e. at temperatures both above and below 
transition. One can see that values of A^(p) are con- 
sistent with zero for all /3-values shown in the figure. 
This was observed on all other lattices employed in 
our study. 

In contrast, the Gribov copy dependence of the 
transverse propagator is rather strong. In Figure [3] 
we show the momentum dependence of At(p) on var- 
ious lattices at /3 = 2.35 (T/T c — 1.1). One can see 
that for fixed physical momentum p the effect of Gri- 
bov copies tends to decrease with increasing volume. 
Such behavior is in agreement with the absence of the 
Gribov copies effects (within Gribov region) in the in- 
finite volume limit, suggested by Zwanziger [30]. On 
the other hand, on given lattice there are always 3 or 
4 minimal values of momentum for which these effects 
are substantial. In particular, At(p) varies between 
0.35 and 0.55 for p = 0, between 0.09 and 0.12 for 
\p\ = Prnin = (2/a) sin (ir/L s ) and between 0.03 and 
0.04 for momentum next after p m %n- Similar observa- 
tions were made at zero temperature as well [25J . 

In Figure U we show the parameter A^ (p) for two 
temperatures, T/T c = 1.1 and T/T c = 2, on lattices 
with approximately equal physical volumes. One can 
see that there is rather weak (if any) dependence of 
the Gribov copy effects on the temperature T. 
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FIG. 2: The momentum dependence of At(p) for three 
temperatures. 
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FIG. 3: The momentum dependence of Ar(p) at j3 = 2.35 
on various lattices. 
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FIG. 4: The momentum dependence of At(p) on two lat- 
tices with approximately equal physical volume. 



of D^p (0) and D T C (0) demonstrate a tendency to de- 
crease; moreover, our data, especially for T = 2T C , 
suggest that (0) and D T C (0) seem to (slowly) con- 
verge in the limit L s — > oo. This convergence is again 
in accordance with a conjecture made by Zwanziger 
in [30| and in accordance with the zero-temperature 
case studied numerically in (22, [25[ . 

On the other hand, our data also suggest that Dt(0) 
is non-zero in the infinite volume limit for both values 
of T, in disagreement with the suggestion made in 
[![. There is no indication for a vanishing transverse 
propagator at zero momentum for increasing volume, 
similarly to the situation for the zero-temperature case 
(22I . |25| . This is in agreement with the refined Gribov- 
Zwanziger formalism [3l|, HH . 

Let us note that one still cannot exclude that there 
are even more efficient gauge fixing methods, superior 
to the one we use, which could make this decreasing 
more drastic. 



The comparison of the results for Ar(p) obtained 
on lattices 4 ■ 32 3 at /3 = 2.512 and 6 x 48 3 at (3 = 
2.635 (both corresponding to the same temperature 
T/T c = 2 and the same physical volume) shows that 
the Gribov copy effects depend weakly on the lattice 
spacing a. 

Let us note that results for At(p) discussed above 
are obtained with the gauge fixing algorithm we have 
chosen, i.e., FSA. The value of At(j>) will be essen- 
tially higher if one uses the OR algorithm to compute 
the fc propagator D T c (p). 

In Figure [5] we show the l/aL s dependence and 
Gribov copy sensitivity of the zero-momentum trans- 
verse propagator Dj<(0) in the deconfinement phase 
(for T = 1.1T C and for T — 2T C ). The difference be- 
tween be values (filled symbols) and fc values (open 
symbols) is rather big, as has been already discussed 
above. However, with increasing size L s the values 



IV. VOLUME AND TEMPERATURE 
DEPENDENCE AND FINITE SPACING 
EFFECTS. 

A. Volume dependence 

As is well-known, the finite volume dependence is 
very strong near the second order phase transition 
point T c . For the transverse propagator Dt(p) this 
dependence can be seen from Figure |6] and for the 
longitudinal propagator Dl (p) from Figure [71 both 
calculated at (3 = 2.3 (slightly above T c ) on various 
lattices. 

Deeper inside in the deconfinement phase the finite 
volume effects are much less pronounced, at least at 
non-zero values of momentum. In Figure |8] we show 
the momentum dependence of the transverse propa- 
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FIG. 5: The l/aL s dependence and Gribov copy sensitiv- 
ity of the transverse propagator Dt(O) at /3 — 2.35 (cir- 
cles) and P — 2.512 (squares). Values of L s are given in 
Table [T] Filled symbols correspond to the be ensemble, 
open symbols to the fc ensemble. The lines are to guide 
the eye. 
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FIG. 8: The momentum and volume dependence of the 
transverse propagators Dt(p) at /3 = 2.35 and /3 — 2.512. 



FIG. 6: The momentum dependence of the transverse 
propagators Dt(p) on various lattices at /3 = 2.3. 



gators Dt{p) on various lattices in the deconfinement 
phase at T = 1.1 T c and T = 2T C . Apart from the 
strong volume dependence at p — 0, we see that finite 
volume effects are rather weak, the sizable effects arc 
seen only for the minimal non-zero momentum p m i n 
on a given lattice. 

The volume dependence of the longitudinal propa- 
gators Dl{p) at T = 1.1T C and at T = 2T C is pre- 
sented in Figure GO Evidently, the volume dependence 
of Dl(p) is even weaker than that of Dx(p), it is rather 
weak even at p = 0. Note that at T = 1.1T C D L (0) 
slowly increases with increasing volume, contrary to 
decreasing of Dt(0). 

For every lattice we observed a well-pronounced 



maximum of the transverse propagator (p) at non- 
zero momentum p$ with |po| ~ 0.4 -j- 0.5 Gev. Deep 
in the deconfinement phase the maximum has been 
found before 

nana in. But we observed it for the 

first time at T = T c and T < T c as one can see in Fig- 
ureQj)] In our recent papers [HLUE] we reported about 
the existence of the maximum of the scalar propaga- 
tor D(p) at non-zero momentum on the symmetric 
lattices (zero-temperature case) when lattice size is 
big enough. Therefore, we conclude that the trans- 
verse propagator Dt{p) has its maximum at po 
for all temperatures, and the behavior of the trans- 
verse propagator in the deep infrared is not consistent 
with the simple pole-type behavior both in confine- 
ment and deconfinement phases. 

It is instructive to compare our results for Dt(p) at 
T = 1.1T C on L s = 48 lattices with respective results 
of Ref. [lj| obtained at this temperature on the lat- 
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FIG. 9: The momentum dependence of the longitudinal 
propagators Dl(p) on various lattices at /3 = 2.35 and 
P = 2.512. 
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FIG. 10: The momentum dependence of the transverse 
propagators Dt{p) near T c on the 4 x 48 3 lattices. 



tices with L4 = 4 and L s = 46 and presented in their 
Figure 1. To make this comparison we made renor- 
malization at /i = 2 GeV as it was made in [l5| . In the 
infrared region we found both qualitative and substan- 
tial quantitative disagreement between our results and 
results of Ref. [15j. In particular, the clear maximum 
which we see in our Figure [5] at T = 1.1T C can not 
be seen from Figure 1 of Ref. . These differences 
between our results and results of Ref. , which are 
just Gribov copies effects, are more essential than dif- 
ferences between our be and fc results discussed in 
Section HTU 

In contrast with Dt{p), the longitudinal propagator 
Dl(p) does not show any trace of maximum at \p\ ^ 
both above and below T c , as one can see in Figure E3 
and Figure [TT] in agreement with results of Refs. [lfj, 
[l5j . This gives an idea that it can be fitted by the 
pole- type behavior (see Section [Vj) . The pole-type 
behavior at high temperature is not surprising since 
at high enough temperature the effective theory is the 
Higgs 3d theory with A4 playing a role of the Higgs 
field. 



B. Temperature dependence 



Thus the transverse gluons in Landau gauge are not 
directly related to confinement [T3, HI[ • 

In contrast, the longitudinal propagator Dl(p) 
demonstrates a drastic jump in its values in the in- 
frared when the critical temperature is crossed from 
above (see Figure ITT]) . The respective renormalization 
constants computed at /1 = 2 GeV differ by a few per- 
cent in this case and this difference has a tendency 
to decrease with increasing /i. The difference in the 
renormalization constants only slightly alters the tem- 
perature dependence of Dl(j>) after renormalization. 
Therefore, D^(p) at small momenta can be considered 
as an order parameter signaling the phase transition. 

Note that the study of the related quantity A^2 = 
{g 2 A 2 E — g 2 A 2 M ) suggests that in the vicinity of T c 
the temperature dependence of Dl(p) can have rather 
nontrivial (non-monotonous) character (34|. Further 
studies at T very close to T c are necessary to clar- 
ify this issue. Let us note also that for the first 
time the fast change of the longitudinal propagator 
near the transition point has been observed in [331 ] (in 
the SU (3) case). This fast change has been recently 
demonstrated in Ref. [HI both in SU(2) and SU(3) 
theories. 



The temperature dependence of the transverse 
propagator Dt{p) near the critical point T c is very 
smooth. Figure [TU] makes comparison of the momen- 
tum dependence of Dt{p) for three temperatures : 
T = 0.9T C , T = T c and T = 1.1T C . Indeed, there 
is no sign of sensitivity to the phase transition. It is 
worthwhile to note that the effect of renormalization 
do not alter this conclusion since the renormalization 
constants computed at /i = 2 GeV differ by less than 
0.5%. 



Deep into the deconfinement phase the decreasing 
of the transverse propagator Dt (p) at p ~ with in- 
creasing temperature (see Figure H]) is in a qualitative 
agreement with dimensional reduction since according 
to dimensional reduction at high temperature Dt(p) 
is to be proportional to (g 2 (T) T)~ 2 . The quantitative 
agreement is not yet expected at temperatures consid- 
ered here. Similarly, the electric propagator Z?l(p) for 
small momenta decreases fast with increasing temper- 
ature, see Figure [S] 



8 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
Ipl; Gev 
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gluon propagators Dl(p) near T c on the 4 x 48 3 lattices. 
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FIG. 12: The momentum dependence of the transverse 
propagator Dt(p) on two different lattices corresponding 
to the same temperature and the same physical volume. 
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FIG. 13: The momentum dependence of the longitudinal 
propagator Dl(p) on two different lattices corresponding 
to the same temperature and the same physical volume. 



smaller spacings) as compared to that employed in 
Ref. [Toj | . This explains the fact that the discretiza- 
tion effects we find are much smaller than in that pa- 
per. 

Contrary to the transverse propagator, the data 
for the longitudinal propagator Dl(j>) show substan- 
tial scaling violations in the infrared. This can be 
seen from Figure Q2] where propagators D l (p) com- 
puted on the same lattices as used in Figure [12] are 
depicted. The renormalization constants computed 
at n = 2 Gev differ only by 4% so the scaling vio- 
lations do not disappear after renormalization. Evi- 
dently, the discretization errors are large at momenta 
smaller than 1 GeV. To reduce the finite cut-off effects 
one should increase L4 or use improved lattice action 
(see, e.g. [13]). 

The results of our study of the scaling behavior at 
zero temperature [IE] suggest that for smaller values of 
/3 (i.e., f3 = 2.35, (3 = 2.3 and (3 = 2.26) discretization 
errors are substantial for both propagators. 



To estimate discretization errors at T = 2T C we 
calculated transverse and longitudinal propagators on 
two different lattices corresponding to the same phys- 
ical 3d volume (ai s ) 3 but with different lattice spac- 
ing. These are lattices L4 = 4, L s = 32, (3 = 2.512 and 
L 4 = 6, L s = 48, (3 = 2.635. In Figure [T2l we show the 
momentum dependence of the transverse propagator 
Dt{p) for these two lattices. One can see good agree- 
ment between results obtained on these lattices for all 
included momenta. This implies that at T = 2T C the 
discretization effects are small even on lattices with 
L4 = 4. We expect that this is true also at higher 
temperatures. 

Let us note that at this temperature we have 
employed larger L4-values and larger f3- values (i.e., 



V. ON THE SCREENING MASSES 

One of the interesting features of the high temper- 
ature phase is the appearance of the infrared mass 
scale parameters : itie ('electric') and hum ('mag- 
netic'). These parameters (or 'screening masses') de- 
fine screening of electric and magnetic fields at large 
distances and, therefore, control the infrared behav- 
ior of Dl(p) and Dx(p)- The electric screening mass 
the has been computed in the leading order of per- 
turbation theory long ago: m? E — ^g 2 T 2 for SU{2) 
gluodynamics. But at the next order the problem of 
the infrared divergencies has been found. On the other 
hand, the magnetic mass uim is entirely nonperturba- 
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tive in nature. Thus a first-principles nonperturbative 
calculations in lattice QCD should play an important 
role in the determination of these quantities. 

As has been already mentioned above, the momen- 
tum dependence of the longitudinal propagator Dl(p) 
in the deep infrared is expected to fit the pole-type 
behavior. Indeed, as an illustration, in Figure [14] we 
show the momentum dependence of the inverse prop- 
agator Dl l [p) at T = 0.9T C and T = 1.1T C . Since the 
volume effects are small enough, at least, at p ^ 0, we 
use data obtained on all lattices listed in Table U with 
exception for p = 0. For this momentum we included 
data for the largest lattice only. 

One can see that at small momenta the dependence 
on p 2 is linear. Thus in the infrared region we have 
used the fitting formula 
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(12) 

The results of the fits are presented in the Table |TTJ 
At T = 2T C we can compare our results with results 
of Ref.[l2| shown in their Figure 3. We find that the 
computed on lattices with £4 = 4 is in good agree- 
ment with respective result of Ref.jl2[. For finer lat- 
tice spacing (L4 = 6) our value is only slightly smaller 
than the value for L4 = 4 (see Table Ql]) indicating 
small scaling deviations, while in Ref. jl2| the value 
obtained on the lattice with L4 = 8 was substantially 
higher. We believe that due to simplicity of the fit- 
ting function p"2")) extracting of from Dl(p) al- 
lows more precise measurement of this quantity than 
its determination from the correlator De{z) used in 
Ref.[ll, . 

In Table [TT] we also show the maximal momenta 
Pmax included into a fit. This value was defined by 
condition that the x 2 /dof value for the fit was smaller 
than 2. One can see that p m ax increases with T. 



p 


T/T c 


niE [Gev] 


m E /T 




2.260 


0.9 


0.41(1) 


1.53(4) 


0.15 


2.300 


1.0 


0.46(1) 


1.54(4) 


0.10 


2.350 


1.1 


0.73(2) 


2.06(6) 


0.30 


2.512 


2.0 


1.21(2) 


2.02(4) 


1.0 


2.635 


2.0 


1.15(3) 


1.92(6) 


1.0 



TABLE II: Values of the screening mass tue obtained from 
fits to eq. (|12p and maximal momenta used in the fit p mnl . 



In contrast, the transverse propagator Dt{p) has a 
form which is not compatible with the simple pole- 
type behavior, so for hum another, different from pole 
mass, definition is necessary. We applied two fitting 
functions to our data for Dr{p)- One of them, Gaus- 
sian function with shifted argument 



FIG. 14: The momentum dependence of the inverse longi- 
tudinal propagator 1/Dl(p) at ft — 2.35 and /3 = 2.26. 



has been used recently in Ref. |25j to fit the T = 
gluon propagator in the infrared region of momenta. 
In ()13[) jtim is a massive parameter, \po\ is momentum 
shift and C a normalization constant. Another fitting 
function is shifted pole propagator of the form 



fp(p) = 



C 



(m*, + (|p|-|Po|) 2 ) 



(14) 



(we keep same notations for fitting parameters). The 
zero momentum was excluded from the fitting range, 
the maximal momentum was determined by require- 
ment that respective X 2 /Ndof was smaller than 1. We 
found that both fits work well in the infrared with 
better performance (larger range) for the fit function 

(Ml. 



p 


T/T c 


mM [Gev] 


m M /T 


Pmax [GeV] 


|po|[Gev] 


2.350 


1.1 


0.56(1)(4) 


1.59(3)(12) 


1.3 


0.40(1) 


2.512 


2.0 


0.78(1)(7) 


1.30(2) (11) 


1.3 


0.51(1) 



TABLE III: Values of the mass parameter mu obtained 
from fits to eq. (|14p and maximal momenta used in the fit 

Pmax - 



(13) 



In Table [TTT] we show fitting parameters for fit func- 
tion eq. qi4[). The second error for vtim is the differ- 
ence from result for fit to eq. (|13[) in which case vtim 
was bigger. The difference in values of |po| for two fits 
was less than 1%. 

We can compare our value for mjn at T/T c = 2 
with result from Ref. [12] presented in Figure 4. In 
that paper result for T/T c = 2 was obtained on 32 2 x 
64 x 8 lattice. They found m M /T = 2.0(3), i.e. much 
higher than our value. Apart from difference in the 
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definition of tum this deviation might be explained by 
the Gribov copy effect, which is much stronger for tum 
than for m^. We need to make computations at higher 
temperatures to make more detailed comparison with 
results of Ref.[H|. 



VI. CONCLUSIONS 

In this work we studied numerically the behavior of 
the Landau gauge longitudinal and transverse gluon 
propagators in pure gauge SU (2) lattice theory in the 
infrared region of momentum values. The special ac- 
cent has been made on the study of the dependence 
of these 'observables' on the choice of Gribov copies. 

The simulations have been performed using the 
standard Wilson action at temperatures from 0.9T C 
up to 2T C on lattices with = 4 and spatial linear 
sizes up to L = 48. For T = 2T C simulations were re- 
peated on lattices with L 4 = 6. For gauge fixing gauge 
orbits enlarged by Z(2) flip operations were consid- 
ered with up to 5 gauge copies in every flip-sector (in 
total, up to 40 gauge copies). The maximization of 
the gauge functional was achieved by the simulated 
annealing method always combined with consecutive 
overrelaxation. 

Our findings can be summarized as follows. 

Similarly to the gluon propagator at T = 0, the 
Gribov copy dependence of the transverse propaga- 
tors Dt{p) is very strong in the infrared, more pre- 
cisely, at a few minimal (for given lattice) momenta. 
At the same time for fixed physical momentum p the 
effect of Gribov copies decreases with increasing vol- 
ume in agreement with [30] . We found no dependence 
of the Gribov copies effects on the temperature or lat- 
tice spacing. 

The Gribov copy dependence of the longitudinal 
propagators Z?l(p) is very weak, at least, at non-zero 
momenta, and is comparable with the statistical er- 
rors (so called 'Gribov noise'). 

We have to emphasize that our conclusions for Gri- 
bov copies effects are relevant for our gauge fixing al- 
gorithm and they might change if another, less effi- 
cient, algorithm is used. 

With increasing size L s the be -values of Dt (0) and 
fc -values of Dt(0) demonstrate a tendency to de- 
crease; moreover, D^f(0) and D S T C (0) seem to (slowly) 
converge in the limit L s — > oo which is in accordance 
with a conjecture made by Zwanziger in [3(| and in ac- 
cordance with the zero-temperature case studied nu- 
merically in [22|, HH|. However, -Dt(O) is non-zero in 
the infinite volume limit, in disagreement with the 
suggestion made in [l|. 

We observed the existence of the maximum of the 
Dt{p) at momenta \p\ ~ 0.4 0.5 Gev not only in 



the deconfinement phase but also for T < T c . Thus 
we confirmed that there is no possibility to explain 
the IR-behavior of the transverse Dt{p) gluon prop- 
agator on the basis of a simple pole-type behavior 
~ l/(p 2 + m 2 ). Instead we fitted this propagator 
to fitting functions eq. flTSl and eq. (fl"4]) with massive 
parameter uim- mjf/T is slowly decreasing with in- 
creasing temperature. To check if this decreasing is 
compatible with gT behavior, as expected for the mag- 
netic screening mass, as well as to compare our results 
for this parameter with results for magnetic screening 
mass, obtained by other authors, we need to repeat 
our computations at higher temperatures. 

For Dl(p) we found good agreement with 'pole- like' 
behavior at small enough momenta p < p max with 
Pmax increasing with T. Our value for tue at T = 
2T C agrees well with result from [12| obtained also on 
lattices with L4 = 4. Again, we need results at higher 
temperatures to compare with other results and with 
the perturbation theory predictions. We shall note 
that our method of computing the in the momentum 
space rather than in the coordinate space gives rise to 
higher precision. 

Away from the transition temperature the longi- 
tudinal propagators Dl{p) demonstrate very weak 
volume dependence. The volume dependence of the 
transverse propagators Dt(p) is strong at p = and 
it is weak at p > 0. 

We found very small scaling violations for Dt{p) at 
T = 2T C comparing results obtained on lattices with 
L4 = 4 and 6. In opposite, for Z?l(p) scaling violations 
in the infrared are substantial. 

We confirmed the observation made in Ref. [l5| 
that the longitudinal propagator Dl(p) in the in- 
frared increases fast when temperature crosses transi- 
tion from above. From results presented in Ref. [l5[ 
for T < T c it is clear that D^(p) has a maximum 
near to T c . This is also in agreement with findings for 



A 2 



(9 A E — g A M ) [3 
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